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Abstract — In this paper, a synthesis method for distributed 
estimation is presented, which is suitable for dealing with large- 
scale interconnected linear systems with disturbance. The main 
feature of the proposed method is that local estimators only 
estimate a reduced set of state variables and their complexity 
does not Increase with the size of the system. Nevertheless, the 
local estimators are able to deal with lack of local detectability. 
Moreover, the estimators guarantee j5^-performance of the 
estimates with respect to model and measurement disturbances. 

I. Introduction 

Estimator design has been an essential part of controller 
design ever since the development of state-space based con¬ 
trollers. A milestone was laid by the Kalman Filter in 1960 

[I] . While in the classical estimator design one estimator 
is used for the entire system, decentralized and distributed 
estimators have gained attention since decentralized control 
became one of the mainstream topics of interest in control 
theory [2], [3]. More results on decenttalized estimation were 
presented in [4], [5], [6], [7]. In decentralized estimation, 
typically, a set of estimators are employed to create estimates 
of local subsystem states with only limited assistance from 
each other. Couplings between those subsystems are treated 
as undesirable disturbances. An important requirement of this 
approach therefore is that the local subsystems are detectable 
from local measurements. 

On the other hand, distributed filtering techniques like 
Distributed Kalman Filtering were presented in [8], [9], [10], 

[II] . In a distributed estimation setup, multiple estimators 
create an estimate of the system’s state, while cooperat¬ 
ing with each other. Situations are not uncommon where 
individual estimators are unable to obtain an estimate of 
the state on its own and cooperation becomes an essential 
prerequisite [12], [13]. At the same time, the progress in 
the area of distributed estimation put forward issues of 
scalability of estimator networks, i.e., there is an interest 
in distributed estimation methods where the dimension of 
the local estimators does not increase with the total size 
of the system. For instance, this is relevant for multi-agent 
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systems, where the agents are not able to perform a self¬ 
measurement, but only receive relative information [14], 
[15]. Direct applications of the existing distributed estimation 
algorithms such as those reported in [8], [9], [12], [13] 
result in the estimators reproducing the entire state of the 
complete network, and therefore, the order of the estimators 
grows with the size of the network. However, for cooperative 
control of multi-agent systems, only local information is 
required for each agent, making the estimate of the whole 
network superfluous. This makes the direct application of the 
mentioned distributed estimation algorithms ineffective with 
respect to the necessary computation power. 

In this paper, our main contribution is the development 
of a new framework that combines the benefits of both 
decentralized and distributed estimation. We develop an 
estimation setup, where local estimators only reproduce a 
desired subset of state variables and their complexity does 
not grow with the total size of the system, in contrast to the 
existing methods for distributed estimation mentioned above. 
Moreover, cooperation between the local estimators will be 
used to deal with possible lack of local detectability, which is 
a major difference to existing results [4], [7]. Therefore, the 
proposed setup is referred to as Cooperative Estimation. In 
particular, we present an J^-based design, which in addition 
provides guaranteed performance with respect to model and 
measurement disturbances. 

This paper significantly extends our results in [16] by 
allowing a more general class of systems. In [16], we 
consider dynamically decoupled agents and take relative 
output measurements. Here, we expand the methodology to 
general large-scale linear systems, where subsystems may 
be physically interconnected. This class of systems is of 
high relevance due to its wide range of applications, such 
as flexible structures [17] and electrical power grids [18]. 

The paper is organized as follows: In Section 2, we present 
some mathematical preliminaries and the system class under 
consideration. Section 3 presents the methodology of coop¬ 
erative estimation design with guaranteed J^:^-performance, 
where the first step is an algorithm for a judicious partition 
of the state space. In order to keep the presentation sim¬ 
ple and focused, we make the simplifying assumption that 
communications between the estimators are perfect. While 
perturbed communication channels may better reflect real- 
life applications, this case can indeed also be considered 
within this estimation setup although the resulting FMI 
conditions become more cumbersome; cf. [19]. Section 4 
illustrates our results with a simulation example and Section 
5 concludes the paper, and gives an overview on future work. 


II. Preliminaries 

Throughout the paper the following notation is used; Let 
A be a quadratic matrix. If A is positive dehnite, it is denoted 
A > 0, and we write A < 0, if A is negative definite. 0 denotes 
a matrix of suitable dimension, with all entries equal 0. 


C. Re-partitioning of the system 

In the next section, we will re-partition the vector x for 
designing local estimators. Associated with the collection 
of outputs 0, for every k = we choose a Oj.- 

dimensional partial state vector 


A. Communication graphs 

In this section we summarize some notation from the graph 
theory. We use directed, unweighted graphs ^ , S') to 

describe the communication topology between the individual 
agents, y = {vi,...,vaj} is the set of vertices, where vt&y 
represents the k-\h agent. S C'Y is the set of edges, 
which model the information flow, i.e. the k-th agent receives 
information from agent j if and only if {vj,Vk) G S. The set 
of vertices that agent k receives information from is called 
the neighborhood of agent k, which is denoted by = 

{j : {vj,Vk) G S^. The outdegree of a vertex k is defined as 
the number of edges in S, which have Vk as their tail. 

B. System model 

We consider a large-scale linear time-invariant system, 
which consists of N interconnected subsystems that are each 
described by the differential equation 

N 

Xk=AkXk-\-Y,^kjXj+BkV, (1) 

./=i 

N 

yk ~ YkXk “t“ ^ ^kjXj “f Bkt (2) 

;=i 

for ^ = 1, ...,A, where Xk G M”* is the state variable, G 
is the output, and v{t),r]i^{t) G-Sf 2 [ 0 ,°°) are .jSf 2 -integrable 
disturbance inputs of subsystem k. The scalar components 
of Xk will be denoted Xk,i. Note that in this system, all 
subsystems are affected by the common disturbance v. This 
assumption does not lead to loss of generality, since it also 
captures the case where the subsystems are affected by 
different disturbances Vk, by simply stacking v<. into one 
vector. 

The global interconnected system can be written as 

x = Ax + Bv, y = Cx-\-rj (3) 

with 
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by using the stacked state and disturbance vectors .x = 
gR" and 77 = [rij, . 

Assumption 1: The global plant (A,C) is observable. 

It is well known that Assumption 1 is a sufficient condition 
in the centralized case. In this paper, this assumption is to 
setup a basic framework under which the state estimation 
problem under consideration is meaningful. 


x« = 


Hk(i) 




(4) 


where ^k{‘) is a selection function that determines, which 
scalar components xy; are included in x^^\ This represents 
a degree of freedom in the design of the estimators and 
all elements of the global state vector x may be chosen 
that are relevant to subsystem k. For instance, x^l may 
contain x^., but it does not have to include all of them, if for 
subsystem k, some parts of its own state are not important. In 
particular, it is required that all xjj which contribute towards 
yii are included in x^^\ As a result, every output can be 
equivalently expressed as 


y, = C«x« + t7,. (5) 


One possible choice of x^^) is the stacked vector including 
Xk and all xj with Ckj ^ 0. In that case. 
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and the rest of the coefficients in ® are defined in a similar 
fashion. 

The selection function is a discrete injective map 

^k-{^,-,Ok]^‘3t', ak<n, ( 6 ) 

where the set = {[kf)\k = \,...,N\i = 1,...,n^} is defined 
as the combination of all appearing indexes of the subsystem 
states and their scalar components x^ ,. For the ease of 
notation, we refer to the elements of the set as A, i.e., 
A = (k,0G^. 

The image of is denoted as c and the 

inverse map is an enumeration of the elements of iS) ^ 

(7) 

which assigns a position in x^^^ to selected components xx 
of the global state vector x. 

In general, partial state vectors x^l may overlap, e.g. x^'^ 
and x*^^^ may contain a common component xx ■ 

For all k = 1,A the global interconnected system 0 can 
now be written as 
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by permutation of the states. 

For every k, the composition of the matrices A bS\ etc., 
is determined by the composition of the partial state variable 
x^^h in turn, the latter is determined by the components of 
the global state x which estimator k seeks to obtain (see 
Section HD- 


















III. Cooperative Estimation Problem and the 
Estimator Design 


The problem considered in this paper is to design a local 
estimator for every subsystem k that creates an estimate for 
the local partial state variable using the local measure¬ 
ments jk described in (|5ll. The vector of local estimates will 
be denoted 



#) = 


.(A:) 


(A:) 

where is the estimate for xx computed at subsystem k. 
The local estimation error vector is defined as 


eW=xW-#) = 


HkW-HkW 


We now formally pose the estimator synthesis problem. 

Problem 1: Determine a collection of estimates xW(f), 
k= such that the following two properties are 

satisfied simultaneously. 

(i) In the absence of model and measurement disturbances 
(i.e., when v = 0, rj =0), the estimation errors decay 
so that eW ^0 exponentially for all k = 

(ii) The estimators (fTTT l provide guaranteed JifL perfor¬ 
mance in the sense that 


<L/ {(^^hf + f\\nkf)dt+Io, 
1 0 


(9) 


for a positive semi-definite weighting matrix 
In (|9]l, /o = is the cost due to the 

observer’s uncertainty about the initial conditions of 
the subsystems. 

This will be achieved by allowing certain agents to com¬ 
municate with each other. 


A. Communication Requirements 

There are two factors, which influence the required com¬ 
munication for the cooperative estimation setup: The first 
one is detectability of In the special case of 

= 0, and being detectable for all A: = 1, ...,W 

no communication is necessary at all, as for every subsystem, 
an estimator can be designed separately. However, these 
assumptions may not hold in a general case. In particular, in 
this paper, we do not require that {A'^^\C^^'^) are detectable 
for all k— 1, ...,W which is a major difference compared to 
existing methods in literature, for instance, [4] and [7]. In 
fact, even all may be undetectable. 

The second factor which influences the required commu¬ 
nication is sparsity of A^^^. Ideally, when the partial state 
is decoupled from the rest of the system, i.e. A^^^ = 0, 
a standard filter can be employed to carry out the 


estimation of x^^^ from y^. However, if x^^^ includes a state 
xx, which is connected to a state x;^* that is not a component 
of x^^\ then the problem becomes more challenging. When 
the connections strength is limited, this can be handled 
using methods like the Small Gain Theorem. Otherwise, 
communication may be used to compensate for coupling 
between x^^^ and x[^^. In this paper, we investigate the latter 
method. 

In order to define the required communication channels, 
we use an assignment function 

C:^r^{0,l,...,iV}, (10) 

with the property that 

if 1^(1) 7 ^ 0. Moreover, i^(A) = 0 only if A ^ 7^*^^ for all k = 
The map ^{■) assigns responsibilities in estimating 
the system’s states to the subsystems and their local estima¬ 
tors. In general, C(A) is not unique and there is a degree of 
freedom in selection of the assignment function. However in 
the case when x^*^) ’s do not overlap, the assignment function 
i^(A) is unique. 

With the definition of the assignment function i^, we can 
introduce the assumption on the communication graph used 
in this paper: 

Assumption 2; If a component xx of x^^^ is physically 
coupled to a state xx*, where A* G then subsystem 

j = C{^*) 7^ 0 can communicate to subsystem k, i.e. {j,k) € 

We denote with Ic^'^ the set of all indexes A* G with 

the property that for all x;^* G Ic^\ there exists a component 
Xx of x^^\ which is coupled to xx*- Assumption 2 reflects 
the point made above, as the more entries A*^*) has, the more 
communication between the subsystems is required. Some 
remarks on the realization of this assumption are in order: 

Remark 1: In the literature (see e.g. [4], [17]), it is often 
assumed that the communication topology simply mimics the 
interconnection topology. Assumption 2 is more general in 
this respect as the selection of the assignment function ^ and 
the definition of the partition x^^^ imply a certain degree of 
freedom. In particular, if the partial state x^*^^ is coupled to 
some state xx* which is not a component of x^^\ there are 
two ways for subsystem k to obtain the required information: 

1) One intuitive option for subsystem k is to receive an 
estimate of xx* through communication with another 
subsystem. The transmitting subsystem may be any 
subsystem, which is assigned to estimate xx* by the 
choice of the assignment function The function 
should therefore be selected so that subsystem j = 
i^(A*) can communicate to subsystem k. 

2) An alternative option is to include xx* into the partial 
state x^^), and therefore, assign the task of estimating xx* 
to subsystem k. This circumvents the communication 
requirement of Assumption 2. This method is meaning¬ 
ful, when the local estimator has enough computational 
power to handle additional coordinates. 







Fig. 1. Example of an admissible estimator structure. The nodes in the 
center represent the subsystems, which receive measurement information 
and estimate the pailial states that are written in the outer circles. The 
black arrows represent the communication links between the subsystems. 


Remark 2: In some applications, the communication 
topology is not a fixed system property, but is a design 
parameter. In that case. Assumption 2 can be easily realized. 

Lemma 1: For all A G ^ there exists a /: € {1, ...,A}, such 
that A G/W. 

Proof: Suppose there exists a A G such that A ^ 
for all k = I, By the definition of the partial states 
and the selection function the column of C which 
corresponds to is 0. Moreover, by the definition of we 
have C,iX) =0 and thus, it follows from Assumption 2 that 
there is no partial state x^^'^ that is coupled to xx. Therefore, 
xx is not observable, which contradicts Assumption 1. 

An example for the interconnection structure which satis¬ 
fies Assumption 2 is shown in Figure [T] This communication 
structure pertains to the numerical example presented later in 
Section IV. As noted for the partition (|4li, the vectors may 
overlap. Therefore, including a consensus term whenever 
overlapping estimators can communicate is able to enhance 
estimation performance of the subsystems and in some cases 
even facilitates feasibility of the design conditions which will 
be introduced later in the paper. 


B. Formal definition of estimators 

The estimator dynamics are now proposed for each sub¬ 
system as 


fW =A«x«+LW(y^-C%W)-f Y. k'"’ 
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with initial condition Xq — 0, where 


Aik) 


J A 


( 11 ) 

is the column 


of aW which corresponds to xx and the unit vector e^-\(^x) 

injects the difference x^^ ~^^x Oi-dimensional space. 

Problem 1 can now be particularized using the estima¬ 
tors (fTTI) . 


Problem 1’: Determine estimator gains Li^\ Ri^i in (fTTIl 
such that properties (i) and (ii) of Problem 1 hold. 

In order to solve this problem, we define the extended 
graph which will be used in the analysis of the in¬ 
terconnection structure between the subsystems. Let every 
subsystem be represented by a cluster of nodes, where 
vertex represents the estimator state x^ . The edges of 
If are now determined by Algorithm [T] An example for the 
extended graph is shown in Figure |2] which applies to the 
numerical example presented in Section IV. 

Algorithm 1 ^ 

Result: Graph If ^ 

Set the edge set S’ = {}. 

for k = do 

for A* G finterconnection partj do 

Add the edge to S, where xx is 

the component of xi^i that is coupled to xx*. 

end 

for j G Mi and A G li^i Cili^i (fusion partj do 
I Add an edge (v;^, ) to S. 

end 

end 


The graph generated by Algorithm 1 graphically displays 
the detailed connection structure of the estimation vectors 
xi^i. The out-degree of vertex in the extended graph is 
denoted by q{k,f). This definition will be used to present 
our main results on the design of the hlter gains, which are 
given in the next section. 


C. Filter gains design 
We define the matrices 
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Fig. 2. Example of the extended graph The subsystems are now 
represented by clusters of vertices, where every vertex represents a single 
estimator coordinate. 
























where e p{k) ^ j^cr^xcr* unknown matri¬ 

ces, pW G 1 ^ 0 * xoj, ^ symmetric, positive definite matrix 
and G M is the i-th diagonal element of and 

a are positive constants which will later play the role 
of design parameters. Furthermore, we define p(A) as the 
diagonal element of P^^W) which corresponds to i.e. the 
^^^^)(A)’th diagonal element. 

Next, for all k— we define the matrices 
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where {A*-^, A^'^ ,-.■}= H 7^^). 

With these definitions, we are ready to present following 
theorem. 

Theorem 1: Consider a group of interconnected LTI sys¬ 
tems ([T]i with local outputs (|2l). Problem 1 admits a solution 
in the form of estimators (fTTT) with 

pik) ^ [p{k)y^Q{k)^ 
p(k) _ (y>{k)ylp{k) ^ 

if the matrices G^^^ and P^^\ k = 1,... ,77 are a solution 
of the following LMIs: 


( 12 ) 
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with {jl.k,j 2 ,k,-;jti„k} = y^k- 



Proof: The estimator error dynamics at 

gW =(A«-LWcW)e^*'^- 
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We use a Lyapunov function 


F(£) = ^ gWTpWg^^ 

v(*)(e(*)) 


where ^re the individual components of y(e). 

The Lie derivative of is 

yW(gW) ^2£WTpW(A('') 
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With the filter gains (fT^ and the LMIs (fOT l it can be 
obtained that 

y W (£) =eWT ('gW _ _ n,) gW 
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Summing up the y(^^s, it holds for V that 
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Integrating both sides of (fT4l) on the interval [0,r], we obtain 

N i-T 

V{s{T)) + y / 

k=i-'^ k=i 

As V{e{T)) >0 and with the zero initial conditions of the 
observer states, it follows that 

y f r((o^\\vf + f\\ri(’^'>\\Adt+Io. 

tiJo tiJo V ) 

Letting L —^ oo, this satisfies Property (ii) of Problem 1. 

Moreover, if = 0 and rjk = 0 for all k= 1, then it 
follows from (fl4l i that 

y(e) < -aV, 

which implies that Property (i) of Problem 1 holds. 

■ 

Note that the choice of a determines the convergence 
speed of the estimators, where a larger a enforces faster 
convergence of the estimates. However, larger values of a 
typically lead to higher filter gain values. 

The salient feature of the resulting cooperative estimators 
(fTTb is that these estimators are local and their complexity 
does not increase with the total size of the network. In 
this sense, the method presented in this paper is scalable 
and guarantees ,^^-type performance. In contrast, a direct 
application of the algorithms developed in [8], [9] and [12], 
[13], to the problem considered here would result in the order 
of the estimators growing with the size of the network. Some 
remarks on the solution of the LMIs (O are in order now. 

Remark 3: As it can be seen from the LMIs (ITsT l. the 
solution to design problem presented here involves solving 
coupled LMIs. When the nature of the application allows 
for these LMis to be solved offline, this can be done in a 
centralized manner. The resulting gain matrices can 

then be deployed to the filters, this will ensure that while 
the estimation algorithm is running, the estimators are fully 
distributed. Alternatively, it was shown in [13], that such 
LMis can be solved with gradient descent type algorithms 
that allow distributed implementation. 

Remark 4: As noted before, the choice of the partial state 
vectors (|4|i is not unique. As a special case, the choice = 
.X for all k= 1,...,A yields local estimators similar to [12], 

[13]. 

IV. Simulation Example 

As an illustrative example, we consider a 9-dimensional 
system, which is composed of 4 subsystems 



t 


(a) Partial state vector xd) and its estimates. 



(b) Partial state vector and its estimates. 



(c) Partial state vector x^^ and its estimates. 



(d) Partial state vector x^) and its estimates. 


Fig. 3. Plots of the local states. Black lines represent the actual partial 
state vectors x^), red lines represent the estimated states x^). 
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( 15 ) 



























with X = [xi,i xi ^2 ^ 2,1 -^ 2,2 21 : 2.3 21 : 3 , 121 : 3,2 21 : 4 , 121 : 4,2 and four 
measurements 
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which are associated with the four subsystems. Thus, the 
partial state vectors (|4|i can be chosen as shown in Figure [T] 
Note that although the system has a dimension of 9 and all 
the coordinates of the state vector are tightly coupled, none 
of the estimators need to handle more than three coordinates. 
In particular, as defined in ([8]l is not detectable. 

This emphasizes the benefits of our algorithm, which can 
deal with situations in which subsystems are individually not 
detectable. 

The graphs representing the interconnection topology are 
shown in Figure [T] and |2] After applying the estimators (fTTT) . 
we obtain simulation results shown in Figure [3 All local 
estimators obtain a correct estimation of their respective 
partial state vector. Moreover, property (ii) of Problem 1 is 
satisfied with performance values of y— 11.1 and m = 8.4. 

V. CONCLUSION 

In this paper, we presented a J^-based approach to coop¬ 
erative state estimation for linear interconnected large-scale 
systems, such as multi-agent systems. In order to achieve 
scalability of the estimation setup, we required the local 
estimators to estimate local states only. We establish an algo¬ 
rithm for interconnecting the local estimators, whereby both 
physical couplings and detectability issues can be handled. 
Moreover, design conditions are presented to guarantee J^- 
performance with respect to both model and measurement 
disturbances. 

Further research will include the distributed calculation 
of the filter gains (O. Moreover, an interesting problem to 
consider is concerned with the application of the cooperative 
estimation algorithm to cooperative and decentralized control 
problems. This may yield interesting results with respect to 
the J^-performance that can still be guaranteed. 
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